Impact of general reionization scenarios on extraction of inflationary parameters 
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Determination of whether the Harrison-Zel'dovich spectrum for primordial scalar perturbations 
is consistent with observations is sensitive to assumptions about the reionization scenario. In light 
of this result, we revisit constraints on inflationary models using more general reionization scenarios. 
While the bounds on the tensor-to-scalar ratio are largely unmodified, when different reionization 
schemes are addressed, hybrid models are back into the inflationary game. In the general reionization 
picture, we reconstruct both the shape and amplitude of the inflaton potential. We find a broader 
spectrum of potential shapes when relaxing the simple reionization restriction. An upper limit of 
10 16 GeV to the amplitude of the potential is found, regardless of the assumptions on the reionization 
history. 
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I. INTRODUCTION 

The inflationary paradigm seems to be the ideal mech- 
anism not only for solving cosmological paradoxes such 
as the observed large-scale smoothness and spatial flat- 
ness of our universe, but also for providing the initial 
seeds for structure formation. The simplest inflationary 
model makes use of a single scalar field (f> (the inflaton), 
which slowly evolves in a very shallow, nearly constant, 
potential V((f>). The dynamics of slow roll gives rise to a 
quasi-de Sitter phase of exponential expansion in the very 
early universe. Since the slope of the spectrum is closely 
related to derivatives of the field potential, slow-roll dy- 
namics predicts that the perturbation spectra should be 
very close to scale invariant, although not exactly so. 
This implies that a quite general inflationary prediction 
is that the power spectra of both scalar, Pr, and ten- 
sor, Pt, fluctuations can be well approximated by power 
laws, i.e., 



Pn(k) oc k 1 



P T {k) oc ¥■ 



(1) 



where the spectral indices n and ut have very mild, if 
any, dependence on the scale k. A scale-invariant scalar 
power spectrum corresponding to the value n — 1 is the 
model proposed by Harrison, Zel'dovich, and Peebles [l[ . 
In other words, from all the considerations above, infla- 
tion predicts n ~ 1, but usually n ^ 10 

A value of the spectral index n slightly different from 
unity would strongly point to the inflationary paradigm 
as the mechanism responsible for providing the initial 



1 For a discussion of slow-roll inflation models with n = 1, see Ref. 
0- 



conditions for structure formation. In addition, in many 
inflationary models the amplitude of gravitational waves 
is proportional to |n — 1|. Confirmation of a deviation 
from a scale-invariant power spectrum would encourage 
the gravitational waves hunters to keep searching for the 
detection of a nonzero tensor amplitude. 

The most recent analysis by the Wilkinson Microwave 
Anisotropy Probe (WMAP) team of their seven-year 
data [3[ rule out the Harrison-Zel'dovich (H-Z) primor- 
dial power spectrum at more than 3cr when ignoring 
tensor modes: n — 0.963 ± 0.012. But this, as well 
as most other previously derived constraints from CMB 
data on cosmological parameters have assumed a "sud- 
den" and complete reionization at a single redshift z r . 
The reionization redshift, z r , is taken to be in the range 
4 < z r < 32, and the cosmological constraints are ob- 
tained after marginalization over z r . The electron ion- 
ization fraction x e {z) is such that for z -C z r x e (z) — 1 
(x e (z) = 1.08 for z < 3 in order to take into account He- 
lium recombination) and x e (z) — 2x 10~ 4 for z > z r , i.e., 
joining the value after primordial recombination with a 
smooth interpolation. 

The process of structure formation that led to gravita- 
tional collapse of objects in which the first stars formed 
are still subject to theoretical and observational uncer- 
tainties. As these first sources began to illuminate their 
local neighborhoods, the HI present in the IGM was 
"reionized." The end of the Dark Ages (the period be- 
tween the end of CMB recombination and the appearance 
of the first stars) remains to be explored and understood. 

There are two main effects on the CMB anisotropies 
produced by the free electrons of the ionized gas: the first 
one washes out the primary anisotropies of the tempera- 
ture autocorrelation (TT) spectrum. The damping of the 
TT signal is quantified by the optical depth parameter 
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r, proportional to the column density of ionized hydro- 
gen. Earlier reionization leads to a the larger suppression 
of the TT acoustic peaks. The second effect produces a 
damping and an additional peak in the polarization au- 
tocorrelation spectrum (EE) The position of this 
new peak in the polarization signal is proportional to 
the square root of the redshift at which the reionization 
occurs, and its amplitude is proportional to the optical 
depth. Since the precise details of reionization processes 
are currently unknown, it is mandatory to explore the im- 
prints of general reionization histories on the CMB spec- 
tra. In the standard, sudden reionization scenario, the 
EE spectrum depends exclusively on the value of Thom- 
son optical depth r. In turn, in extended reionization 
schemes, the precise history of how the universe became 
ionized affects the large-scale EE power spectrum in a 
crucial way [5j, and the power is transferred from larger 
to smaller scales when considering that reionization pro- 
cesses could take place in a non-negligible redshift (time) 
interval. 

The major goal of this paper is to study how current 
constraints on the scalar spectral index n and the tensor- 
to-scalar ratio r are modified if the standard ("sudden") 
reionization assumption is relaxed. In a precursor study 
we demonstrated that in a general reionization scenario 
the Harrison-Zel'dovich spectrum (n = 1) is perfectly 
consistent with observations |6|. In this study we shall 
also include information from tensors modes, showing 
that inflationary models that are ruled out in the sudden 
reionization scheme are allowed in more general reioniza- 
tion scenarios. We also reconstruct both the shape and 
the amplitude of the inflationary potential V(<f>) allowed 
by current data in both sudden and general reionization 
schemes. 

The paper is organized as follows. Section [TT] sum- 
marizes the results of inflationary theory relevant for our 
considerations. A possible classification of different mod- 
els of inflation is presented in Sec. IIII1 The analysis 
method used here to derive the cosmological constraints 
is described in Sec. US Section [V] gives the resulting 
constraints on cosmological parameters and their impli- 
cations for inflationary models. The inflationary poten- 
tial reconstruction method and the results are presented 
in Sec. EU We conclude in Sec. IVHl 
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potential V(<f>) with equations of motion 



INFLATION AND THE HAMILTON- J ACOBI 
FORMALISM 



In this section we briefly review the dynamics of a 
scalar field in a cosmological background. We assume 
a flat Friedmann-Robertson- Walker (FRW) metric: 



H 2 = 



3H<p + 
Sn 



3r7ipj 



V(<f>) 
1 
2 







-0 2 + y(0) 



(3) 
(4) 



where H = a/ a is the Hubble parameter, dots and primes 
denote derivatives with respect to cosmological time and 
to the scalar field respectively, and mpi is the Planck 
mass. From the definition of the Hubble parameter, it 
follows that 



a(t) oc e 



-JV 



exp 



H(t)dt 



where the number of e-folds N is simply 



N = 



H(t)dt 



(5) 



(6) 



where t e refers to the end on inflation. The integration 
cxtrcma are chosen in such a way that N = coincides 
with the end of inflation. 

A very powerful way of describing the inflationary dy- 
namics is given by the Hamilton- Jacobi formulation of 
inflation. The basic idea is to consider the scalar field cf) 
itself to be the time variable; this can be done as long as 
it varies monotonically with time. Then, expressing the 
Hubble parameter as a function of the field, H = H(cf>), 
the equations of motion become 



m P1 



32tt' 



V(4>). 



(7) 
(8) 



'pi 



The second of these equations is called the Hamilton- 
Jacobi equation. Inflation takes place while the field is 
slowly rolling towards a minimum of the potential, and 
the field energy density is dominated by its potential en- 
ergy. More quantitatively, the slow-roll approximation 



holds in the limit in which <j) <C 3H<f> and 
that Eqs. ([3]) and (gj) become 



< V, so 



H 2 



3H 
8tt 
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v{<t>) 



(9) 



The validity of the slow-roll approximation is quite nat- 
ural because the slope of the inflaton potential must be 
sufficiently shallow to drive inflation. For single-field in- 
flation during the slow-roll phase, the kinetic energy of 
the field is negligible and the potential is nearly constant: 



ds 2 = dt 2 - a 2 (t)[dr 2 + r 2 dn 2 



(2) 



= V(4>) + ^ ~ V(<j>) 



const. 



(10) 



where a(t) is the cosmological scale factor. In an FRW 
background, a scalar field <f> evolves under the action of 



From Eq. ©, we can see that this gives rise to a (quasi-) 
de Sitter phase with H almost constant. The amplitude 
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of the potential must be sufficiently large to dominate 
the energy density of the Universe at that epoch. 

The slow-roll approximation is consistent if both the 
slope and the curvature of the potential are small (in 
units of the Planck mass) when compared to the potential 
itself: V, V" -C V, or equivalently if the so-called slow- 
roll parameters e and 77 are much smaller than unity. The 
slow-roll parameters are defined as 
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(11) 



When V', V" <C V, both e and r\ can be expressed in 
terms of the potential and its derivatives as 
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(12) 



The consistency of the slow-roll condition thus implies 
that e, \r/\ <C 1. Using the definition of e, the Hamilton- 
Jacobi equation can be rewritten in the useful form 



H 2 (<f>) 



1 



1 



8tt 
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V(<f>) 



(13) 



Inflation also provides a natural mechanism to gener- 
ate the inhomogeneities presently observed in the Uni- 
verse. During inflation, quantum fluctuations, inevitably 
present at small scales, are quickly redshifted to scales 
much larger than the horizon size and then frozen in as 
perturbations to the background metric. The perturba- 
tions created during inflation can be of two types: scalar 
(or curvature) perturbations, which couple to the matter 
stress-energy tensor, and tensor perturbations (gravita- 
tional waves), which do not couple to matter. The power 
spectrum of scalar perturbations (quantified as pertur- 
bations in the Ricci scalar TV) is described by 



(*) = 


H 1 " 







and its spectral index n reads 

din fc 



(14) 



(15) 



The power spectrum of tensor fluctuation modes is given 
by 



P^ 2 (k) = 



4 H \ 



(16) 



again evaluated when the mode k crosses the horizon. 
The ratio of the tensor-to-scalar perturbation is defined 

as 



Pr 
Ptz 



(17) 



and, as in the scalar power spectrum case, one can write 
Pt oc k riT . The two spectral indices expressed in terms 
of the slow-roll parameters are 



n 
tit 



l-4e + 2i] 
-2e , 



and the tensor-to-scalar ratio r is 



r = 16e 



(18) 
(19) 



(20) 



The relations above are valid at first-order approximation 
in the slow-roll parameters. Therefore, if primordial per- 
turbations originated from the dynamics of a slow-rolling 
scalar field, the spectrum should not be exactly scale in- 
variant. In fact, since the slow-roll parameters e and 77 are 
small, but not vanishing (in other words, since the poten- 
tial is very close to flat but not exactly flat), we expect 
that n ~ 1 but nevertheless n 7^ 1. A scale-invariant 
power spectrum corresponds to the value n = 1 is the 
aforementioned model proposed by Harrison, Zel'dovich, 
and Peebles [l|. Given the fact that Pr cx fc n_1 , the 
spectral index can be thought as a measure of the depar- 
ture of the spectrum of the scalar perturbations from an 
exactly scale-invariant power spectrum. 



III. ZOOLOGY OF INFLATIONARY MODELS 

In this section we follow the classification of Kinney 
et al. 7}. At lowest order in the slow-roll approxima- 
tion the relevant parameters to distinguish among infla- 
tionary models are n and r Q. The different classes of 
models are characterized by the relation between these 
two parameters, or equivalently, by the relation between 
e and 77. At lowest order in the slow-roll approximation 
we can divide the inflationary models into three general 
types: large-field, small-field and hybrid. The boundary 
between large-field and small-field models is represented 
by the so called linear models. 

• Large-field models are characterized by — e < 
-q < e. Popular examples of large-field models 
are V{(f>) = A 4 (</>//i) p and exponential potentials, 
^) = A 4 e X p(^/V). 

• Small-field models are characterized by rj < — e. 
They result from a generic potential of the form 
V{(j>) = A 4 [l — (<f)/(j,) p ], which can be understood 
as the lowest-order Taylor expansion of an arbitrary 
potential about the origin. 

• Hybrid models are characterized by < e < r\. 
A generic hybrid potential is of the form V(4>) — 

• Linear models are on the boundary between large- 
field and small-field, and they are characterized for 
this reason by rj = —e. The generic linear potential 
is of the form: V(<p) oc 4>. 
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With the above classification we can cover the entire n-r 
plane and derive constraints on the inflationary models 
directly from the constraints on the n-r plane that arise 
from cosmological observations; see Sec. [V] 



IV. ANALYSIS METHOD 

We adopt two different methods for parametrization of 
the reionization history. The first method, developed in 
Ref . [5| , is based on principal components that provide a 
complete basis for describing the effects of reionization on 
large-scale -E-mode polarization. Following Ref. one 
can parametrize the reionization history as a free func- 
tion of redshift by decomposing x e [z) into its principal 
components: 



x e (z) = x{(z) + ^ m n S n( z )> 



(21) 



where the principal components, S^z), are the eigen- 
functions of the Fisher matrix that describes the depen- 
dence of the polarization spectra on the electron ioniza- 
tion fraction x e (z), are the amplitudes of the princi- 
pal components for a particular reionization history, and 
x{(z) is the WMAP fiducial model at which the Fisher 
matrix is computed and from which the principal compo- 
nents are obtained. In what follows we use the publicly 
available S^z) functions and vary the amplitudes for 
fi = 1, 5 for the first five eigenfunctions. Hereafter we 
refer to this method as the MH (Mortonson-Hu) case. 

In a second approach to a general reionization prescrip- 
tion we employ a different parametrization, sampling the 
evolution of the ionization fraction x e as a function of red- 
shift z at seven points (z = 9, 12, 15, 18, 21, 24, and 27), 
and interpolating the value of x e (z) between them with 
a cubic spline. For 30 < z we fix x e = 2 x 10~ 4 as the 
value of x e expected before reionization (and after pri- 
mordial recombination), while x e = 1 for 3 < z < 6 and 
x e = 1.08 for z < 3 in order to be in agreement with both 
Helium ionization and Gunn-Pctcrson test observations. 
This approach is very similar to the one used in Ref. [t| , 
and we will refer to it as the LWB (Lewis- Weller-Battye) 
case. 

We then modified the Boltzmann CAMB code [lol ]. 
incorporating the two generalized reionization scenarios 
and extracted cosmological parameters from current data 
using a Monte Carlo Markov Chain (MCMC) analysis 
based on the publicly available MCMC package cosmomc 

We consider here a flat ACDM universe described by 
a set of cosmological parameters 



{uj bl uj Ci @ s ,n, log[10 10 ,4 s ], r, n run }, 



(22) 



where ojf, = ilbh 2 and oj c = Q c h 2 are the physical 
baryon and cold dark matter densities relative to the 
critical density, S is the ratio between the sound hori- 
zon and the angular diameter distance at decoupling, A a 



is the amplitude of the primordial spectrum, n is the 
scalar spectral index, r is the tensor-to-scalar ratio, and 
n run = dn/dhik is the running of the scalar spectral 
index: 



A^(fc) = A^(fco) 



n(ko) — l+rr ln(fc/ ko)dn/ din k 



(23) 



Here, ko = 0.002 Mpc -1 is the pivot scale. 

The extra parameters needed to describe reionization 
are the five amplitudes of the eigenfunctions for the MH 
case, or the seven amplitudes in the seven bins for the 
LWB case, and one single common parameter, the optical 
depth, t, for the sudden reionization case. 

Our basic data set is the seven-year WMAP data 
Q (temperature, polarization, and tensor modes) with 
the routine for computing the likelihood supplied by the 
WMAP team. Together with the WMAP data, we also 
augment the WMAP7 data with the CMB data sets from 
BOOMERanG 0, QUAD 13], ACBAR [3, and BI- 
CEP [15j. For all these experiments we marginalize over 
a possible contamination from the Sunyaev-Zel'dovich 
component, rescaling the WMAP template at the cor- 
responding experimental frequencies. We therefore con- 
sider two cases: we first analyze the WMAP data alone, 
referring to it as to the "WMAP7" case, and we then 
include the remaining CMB experiments ("CMB- ALL"). 



V. RESULTS 

Table U summarizes the main results of the analysis for 
different cosmological data sets, showing the constraints 
on n and r for the MH, LWB, and sudden reionization 
schemes. When the sudden reionization assumption is re- 
laxed, the mean values of n and r tend to shift to higher 
values. The shift in n was already noted in the previous 
paper [f| . The importance of this shift is that in a gen- 
eral reionization scheme the H-Z spectrum is perfectly 
consistent. Notice, however, that the presence of tensors 
and/or a running spectral index in the analysis allows for 
a H-Z spectrum even in the sudden reionization scheme 
(at a confidence level (c.l.) corresponding to 2a). Nev- 
ertheless, in the case of the MH reionization scenario, 
without running of the index, the best fit for the scalar 
spectral index is already higher than one at 68% c.l. In a 
general reionization scenario the allowed values of r also 
shift to higher values. When additional data from other 
CMB probes are added to the WMAP7 data, the con- 
straints on n and r are shifted back toward lower values. 
In summary, in the MH reionization case ignoring run- 
ning of the spectral index, using WMAP7 data the H-Z 
spectrum (n = 1) is very close to the best fit value, and 
inside the 68% c.l. for the case CMB- ALL. 

The values for the tensor-to-scalar ratio r and the run- 
ning of the spectral index n run at 95% c.l. are slightly 
higher considering a general reionization scenario. How- 
ever, their 68% c.l. constraints barely change when the 
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WMAP7 


CMB-ALL 


Planck 




Sudden 


MH 


LWB 


Sudden 


MH 


LWB 


Sudden 








no 


running of sea 


ar spectral index 
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^Do/o C.I.J 


0.987 ± 0.020 


1.001 ±0.027 


0.992 ± 0.021 


0.974 ± 0.016 


0.985 ± 0.020 


0.977 ±0.017 


0.960 ± 0.004 


n 


[\JD/Q C.I.J 


n < 1.031 


n < 1.067 


n < 1.039 


n < 1.007 


n < 1.026 


n < 1.012 


n < 0.968 


V 


L.I.J 


0.142 ±0.116 


0.141 ±0.119 


0.149 ± 0.115 


0.095 ± 0.079 


0.101 ±0.085 


0.103 ±0.087 


0.053 ± 0.022 


r 


{i)O/0 C.I.J 


r < 0.373 


r < 0.376 


r < 0.371 


r < 0.251 


r < 0.266 


r < 0.275 


r < 0.093 








running of scala 


r spectral index 








n 


(68% c.l.) 


1.067 ± 0.062 


1.080 ±0.065 




1.094 ± 0.052 


1.106 ±0.054 




— 




(QK.% r 1 1 


n < 1.192 


n < 1.207 




n < 1.197 


n < 1.222 


— 




r 


(68% c.l.) 


0.191 ±0.154 


0.200 ±0.158 




0.181 ±0.141 


0.185 ±0.140 




— 


r 


(95% c.l.) 


r < 0.497 


r < 0.515 




r < 0.451 


r < 0.445 






n 


run (68% c.l.) 


-0.040 ± 0.029 


-0.036 ±0.031 




-0.056 ± 0.021 


-0.058 ±0.022 







TABLE I. Constraints for different data sets on n, r, and r Iun in different reionization scenarios with and without the running 
of the scalar spectral index n. 




Hybrid 
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n n 
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FIG. 1. Two-dimensional contour plots at the 68% and 95% confidence levels without running of the scalar spectral index for 
the WMAP7 data (left figure) and the CMB-ALL data set (right figure) . Shaded contours correspond to the sudden reionization 
approximation, while open contours model reionization as MH. The dark solid (dashed) lines refer to large-field models with 
p — 2 (p = 4). The lighter cross (dashed) curves depict small- field models with p — 2 (p = 4). The solid horizontal line that 
basically coincides with the x axis depicts hybrid models with p — 2 (the p — 4 case basically overlaps the p — 2 case). The 
filled circles (squares) denote the points in the parameter space for which the number of e-folds iV is equal to 60 (50). 



reionization history is modified, as expected, due to the 
large uncertainties on r and n run . 

The shift induced on allowed values of inflationary 
parameters n and r by different assumptions for the 
reionization history is important for the subsequent con- 
straints on inflationary models. To study this, we have 



reconstructed the relation between n and r in the differ- 
ent classes of models described in the previous section, 
and we have plotted these relations in the n-r plane, to- 
gether with the cosmological constraints. 

Figure Q] depicts the 68% and 95% c.l. allowed con- 
tours by the WMAP7 data and the CMB-ALL data sets 
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without running of the scalar spectral index for different 
assumptions of the reionization history. The indicated 
contours denote the allowed regions when tensor modes 
are included in the analysis, and when the reionization is 
assumed to be sudden and when using the MH procedure 
(see the figure caption for details). 

Figure [2] is the same as Fig. Q] but now allowing for a 
running of the scalar spectral index. 

Following Ref. Q we can easily develop the different 
expressions concerning the n-r parameter space. For in- 
stance, for large-field models, with a polynomial potential 
V oc <$P , the relation among these parameters is 



r / 2 
1 - - I 1 + - 

P 



(24) 



The dark lines in Figs. Q] and [2] refer to this relation for 
quadratic (p = 2) and quartic (p — 4) potentials. It 
is straightforward to relate r with N (the number of e- 
foldings before the end of inflation) : 



r = Ap/N 



(25) 



which allows us to draw points with N — 50 (squares) 
and N = 60 (circles) in Figs. Q] and El 

Similarly, we can relate n to r, and both of them in 
terms of N for both small-field and hybrid models. For 
small-field models, the generic potential we are using is 
of the form V(cf>) = A 4 [I - (0/^) p ]. Typically in these 
models the slow-roll parameter e (and hence, r) is close 
to zero. The spectral index can be written as 
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It is straightforward to see that for p = 2 
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(26) 



(27) 




0.05 



0.975 



FIG. 3. The 68% and 95% c.l. constraints forecast on the n 
vs. r plane from Planck mock data for n = 0.96 and r = 0.05 
and sudden reionization (dark contour), and (wrongly) fitted 
assuming MH reionization (light contour). 



while for p = A we have 



I - - 



TV 



I6/x 4 



1/3 



3 
N 



(28) 



Figs. Q] and [5] also contain the small-field model case, 
depicted by cross for p = 2 and by the indicated dashed 
curve for p = 4 (assuming /i ps mp{). 

For hybrid models, the potential chosen is V(4>) — 
A 4 [I — a(mpi/<fi) p ], based on potentials generated in dy- 
namical SUSY breaking models [16] . As in small-field 
models, the tensor-to-scalar ratio r is negligible. The 
expression for n is given by 

n«l + 2 7 fetli (29) 

where A^ot is the total number of e-foldings (chosen to 
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be 100 in this example). Notice that Eq. (|29|) indicates 
that the power spectrum in these sort of models is blue 
(n > 1). Indeed, in the sudden reionization scenario 
with negligible running of the spectral index these hybrid 
models are highly disfavored; in more general reionization 
schemes such models are allowed by WMAP7 data; see 
Fig. [TJ When more CMB data sets are included in the 
analysis, hybrid inflation models with a blue tilt are again 
disfavored at 95% c.L, even in the more general reioniza- 
tion scenarios considered here; see the CMB- ALL part 
of Fig. [TJ When a running scalar spectral index is al- 
lowed, hybrid models are perfectly compatible with data, 
regardless of the assumptions about the reionization pro- 
cesses; see Fig. [2] 

The LWB reionization scheme leads to very similar 
constraints to those of MH parametrization on the n — r 
plane (albeit slightly closer to the sudden reionization 
case). Indeed, n is constrained to be red at the 68% c.L 
in the CMB- ALL case, but the H-Z model is still consis- 
tent with data within two standard deviations. 

We also forecast future constraints from the Planck 
experiment with the specifications of Ref. 17], assuming 
that n = 0.96 and r = 0.05 and sudden reionization. 
If the data is (wrongly) fitted assuming a more general 
reionization scenario (MH reionization, for instance), the 
constraints that one would obtain on the n — r plane are 
shown in Fig. [3] Notice that Planck will be able to tell 
n 7^ 1 at a very high confidence level even if the precise 
details of the reionization processes are unknown. Planck 
data will also be sensitive to the tensor-to-scalar ratio at 
the 95% c.L for r > 0.05. 



by the following set of equations [18j, [22l |2E 

2c), 



de 



da n 
— - = -5ecr - 12e 2 + 2A 
dN 



(2) 
H » 



dN 



-2)e 



M+i) 



(31) 
(32) 

(33) 



Given a solution to these equations, the observable 
quantities, i.e., the scalar spectral index n, its running 
n run , and the tensor-to-scalar ratio r, can be evaluated. 
To second order in slow roll, these are given by [2l|, [24j : 

r = 16e[l-C(o + 2e)], 

n-l = cr - (5- 3C)e 2 - 7(3 - 5C)ae 



1 



dn 
dN 



1 



(34) 
0.577 is 



where C = 4(ln2 + 7) - 5 = 0.08145 and 7 
Euler's constant. 

The solution to Eqs. (|31ti33f) also allows one to recon- 
struct the form of the potential V{(f>) (II \^2T\ . In fact, 
from the Hamilton- J acobi equation (see Eq. (fTT]) ) 



v{4>) 



/3m 2 



pi 



V 8vr 



H 2 W 



1 



1 



<<f>) 



(35) 



Once e(N) is known from the solution of Eqs. (|3HI33I) . 
H(N) can be determined from 



VI. MONTE CARLO RECONSTRUCTION OF 
THE INFLATIONARY POTENTIAL 



In this section we describe the technique known as 
Monte Carlo reconstruction, a stochastic method for in- 
verting observational constraints to obtain an ensemble 
of inflationary potentials compatible with observations. 
The method is described in more detail in Refs. [i~8T - [2Cjl | . 

The slow roll parameters e and 77 already have been 
defined in Sec. flJJ see Eq. (fTl) , These two parameters 
are related to the observables n and r by the formulae 
given in Sec. [IIJ valid to first order in the slow-roll ap- 
proximation. We find it convenient to use the parameter 
(7 = 2,7] — 4c in place of 77; the advantage is that to first 
order in slow roll, a ~ n — 1. 

Taking higher derivatives of H, one can construct an 
infinite hierarchy of slow- roll parameters [2l[ : 



w = fgpiV ( g, )^~ 1 d(i+1)H 

4tt ) H e dej)^ 1 ) 



(30) 



1 dH 
HdN 



(36) 



The solution to the above equation allows one then to 
obtain V(N) up to a normalization constant; this is fixed 
by the normalization of the Hubble parameter that enters 
the above equation as a integration constant. We will 
return on this later. 

Finally, in order to obtain 4>(N), we note that Eq. 
and dN/dt = —H together imply that 



'pi 



H' 



TOpi 



dN 



4tt H 2^ 



(37) 



The evolution of the slow-roll parameters is described 



where it should be implicitly understood that there is a 
sign ambiguity in the last equality since it should have the 
same sign as H'(cf>). Since we do not know in advance the 
sign of H', we should consider the two cases separately. 
However, it is easy to see that they are related by the 
transformation </)—>—</>. For this reason, in the following 
we just consider the case dip/dN < (i.e, the value of <j> 
grows as inflation goes on). Once <fi(N) is obtained from 
the solution to the above equation, it can be inverted 
to obtain N(4>) and finally V{4>). However, the value 
of 4>(N) can be known only up to an additive integration 



8 





WMAP7 


CMB-ALL 




Sudden MH 


Sudden 


MH 


n (99% c.l.) 


0.913 < n < 1.221 0.918 < n < 1.234 


0.968 < n < 1.216 


0.984 < n < 1.260 


r (99% c.l.) 


r < 0.665 r < 0.688 


r < 0.551 


r < 0.619 


n run (99% c.l.) 


-0.098 < n run < +0.037 -0.098 < n run < +0.050 


-0.099 < n run < -0.001 


-0.117 < n run < -0.001 



TABLE II. Observational limits (at 99% confidence level) used to constrain the inflationary potential for the four cases considered 
in reconstruction. 



constant. We fix the latter so that (j> = at the beginning 
of inflation. 

In order to calculate an ensemble of potentials that are 
compatible with observations, we proceed in the following 
way: 

1. Choose random initial values for the inflationary 
parameters in the following ranges: 



N 
e 

(7 

(2) 
H 

(3) 



[40, 70] 
[0, 0.8] 
[-0.5, 0.5] 

[-0.05, 0.05] 

[-0.005, 0.005] 



>(6) _ 



truncated at M = 6. 



2. 



Evolve forward in time (dN < 0) until either (a) 
inflation ends (e > 1), or (b ) the evolution reaches 
a late-time fixed point (e = A)+ =0, a = const). 

In case (a), evolve N e- folds backwards in time from 
the end of inflation and calculate the observables 
n — 1, r, and the running n run at that time; in 
case (b), calculate the observables at the time the 
evolution reaches the fixed point. 



4. Repeat the above procedure Nmc times. 

5. Choose a window of acceptable values for the ob- 
servables n — 1, r, and the running n run , and then 
extract from the Nmc models those that satisfy the 
observational constraints. 



the previous section for the "WMAP7" and "CMB-ALL" 
datasets in the two cases of sudden and MH reionization 
scenarios. For convenience, we report four sets of con- 
straints (at the 95% and 99% confidence levels) in Table 
HT1 We show a sample of 300 reconstructed potentials 
in Fig. [4] We have rescaled all the potentials so that 
V(4> — 0) = 1 and < cj) < 1, so that the figure only 
contains information concerning the shape of the infla- 
tionary potential. We see that the WMAP7 data alone 
do not really constrain the shape of the inflationary po- 
tential, even when more general models of reionization 
are considered. However, the situation changes dramat- 
ically when other CMB experiments are included, as it 
can be seen from the right panels of Fig. SJ Presumably 
this is due to the fact that when these experiments are 
included, models with n run = are placed outside the 
observationally allowed region. In this case, the reion- 
ization model makes some difference. In the case of sud- 
den reionization the possible shapes of the potential are 
severely restricted; in particular, it seems that potentials 
with V"{4>) < are not allowed. When a more general 
model of reionization is used, however, more shapes are 
allowed. This is probably an result of the fact that the 
constraints on the scalar spectral index n are weakened. 

Other than the shape of the potential, it is also im- 
portant to constrain its amplitude. The reconstruction 
procedure described above does not yield the amplitude 
of the potential; this has to be fixed from some obser- 
vational input, like the normalization of the Hubble pa- 
rameter. We choose to normalize the Hubble parameter 
through the condition on the density contrast: 



Sp^l H 1 
p 2ir mpi ^/e 



10- 



(38) 



6. Reconstruct the potential for these models, follow- 
ing the procedure described above. 

We have implemented the procedure described above 
with Nmc — 5 x 10 5 . We consider four sets of obser- 
vational constraints, corresponding to those obtained in 



Once this is done, we find that in all the four cases illus- 
trated in Table UH V(4>) < 10 _11 mp 1 . This correspond 
to an upper limit to the energy scale of inflation of about 
10 16 GeV. 



9 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

</> 4> 



FIG. 4. Sample of three hundred reconstructed potentials for the four sets of observational constraints discussed in the text. 
The potentials are calculated at the 99% c.l. All the potentials have been rescaled so that V(cj> = 0) = 1, and < (j> < 1. 
Upper left panel: WMAP7 with sudden reionization; upper right panel: CMB-ALL with sudden reionization; lower left panel: 
WMAP7 with MH reionization; lower right panel: CMB-ALL with MH reionization. 



VII. CONCLUSIONS 



Details of the reionization processes in the late universe 
are not very well known. In the absence of a precise, 
full-redshift evolution description of the ionization frac- 
tion during the reionization period, a simple parametriza- 
tion with a single parameter z r has become the standard 
reionization scheme in numerical analyses. More general 
reionization schemes have been shown to allow values of 
the scalar spectral index consistent with a scale- invariant 
power spectrum. In this paper we deduce information 
about tensor modes, and explore how the inflation con- 
straints are modified when the standard reionization as- 
sumption is relaxed. The tensor-to-scalar ratio bounds 
are largely unmodified under more general reionization 
scenarios. Therefore, present (future) primordial grav- 
itational wave searches are (will be) unaffected by the 
precise details of reionization processes. In the absence 
of a running spectral index, hybrid models, ruled out in 
the standard reionization scheme, are still allowed at the 
95% c.l. by WMAP7 data. The constraints on other in- 
flationary models, such as large-field or small-field mod- 
els, do not change. Future Planck data will be able to 
measure the scalar spectral index n with unprecedented 



precision and be sensitive to tensor modes if r > 0.05 
at the 95% c.l. We also show the impact of different 
reionization histories on the reconstruction of the infla- 
ton potential. The variety of the reconstructed shapes is 
larger in extended reionization scenarios than compared 
to the standard one. Namely, models with V"(<j>) < are 
allowed in general reionization scenarios. However, the 
constraints on the amplitude of the potential remain un- 
changed: we find an upper bound for the latter of about 
10 16 GeV, independent of the reionization process details. 
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